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ABSTRACT 


The  energy  distributions  of  electrons  of  about  54,  75,  and  92 

MeV  were  measured  before  and  after  passing  through  lead  absorbers  of 

2 
thicknesses  ranging  from  0.706  to  2.825  gm/cm  ,  and  through  copper 

2 
and  aluminum  absorbers  up  to  5.726  gm/cm  .   The  electrons  were  acceler- 
ated by  the  LINAC  of  the  Naval  Postgraduate  School.   A  measurement  to 
determine  the  optimum  location  of  the  absorber  indicated  that  the 
separation  between  the  scattering  foil  and  absorber  should  not  exceed 
6  cm,  if  geometric  difficulties  are  to  be  avoided.   The  most  probable 
energy  losses  agree  with  the  theory  of  Blunck  and  Westphal  for  all 

three  materials  and  all  thicknesses.   The  half-widths  of  the  distri- 

2 
butions  agree  with  theory  up  to  thicknesses  of  about  2  gm/cm  .   For 

greater  thicknesses  the  experimental  half-widths  of  copper  absorbers 

appear  to  agree  with  theory;  the  experimental  half-widths  of  aluminum 

are  greater  than  theoretical  predictions;  and  those  of  lead  are  smaller 

than  theory  predicts. 
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I.   INTRODUCTION 

The  energy  distribution  of  an  initially  monoenerget ic  electron 
beam  which  has  passed  through  a  layer  of  matter  has  been  calculated  by 
Blunck  and  Westphal  [l].   This  calculation  evolved  from  the  earlier 
theoretical  work  of  Landau  [2*]»  Eygers[3],  Bethe  and  Heitler  [4],  and 
Blunck  and  Leisegang[5") .   It  essentially  accounts  for  the  energy  loss 
distribution  due  both  to  ionization  and  excitation  and  to  radiation. 
The  theory  is  based  on  the  assumption  that  the  absorbing  material  is 
sufficiently  thin  so  that  the  energy  loss  is  small  compared  to  the 
initial  energy  of  the  electrons.   The  development  of  the  theory  is 
contained  in  Section  II. 

A  number  of  experiments  to  measure  the  energy  loss  distribution 
of  high  energy  electrons  in  thin  layers  of  various  materials  have  been 
performed,  with  incident  energies  in  the  range  10  MeV  to  150  MeV.   Of 
particular  interest  are  the  works  of  Breuer  on  aluminum  at  the 
Darmstadt  linear  accelerator  [6];  Bumiller,  Buskirk,  Dyer,  and  Miller 
on  aluminum  at  the  LINAC  of  the  Naval  Postgraduate  School  [12,16], 
and  Goodwin  on  copper,  also  at  the  Naval  Postgraduate  School  [l3]. 
The  experimental  works  listed  above  are  generally  consistent  in  their 
results . 

In  each  case  agreement  was  found  between  experiment  and  theory  for 

2 
thin  absorbers  (<  2.0  gm/cm  ).   However,  for  thicker  layers  of  material 

the  theory  fails  to  predict  the  energy  distribution  experimentally 

realized.   It  is  not  clear  how  rapidly  the  theory  breaks  down  as  the 

energy  loss  becomes  significantly  large  compared  to  incident  energy. 


In  this  thesis  data  are  presented  for  the  energy  loss  distribution  in 

2 
lead  for  absorber  thicknesses  from  0.706  to  2.825  gm/cm  .   Additional 

data  were  obtained  for  thick  layers  of  aluminum  and  copper  to  supple- 
ment the  experimental  results  obtained  by  Bumiller  et.al.  [15]  and 
Goodwin  [13],   The  energy  losses  ranged  from  \%   to  about  \1°L   of  the 
incident  beam  energy.   The  energy  loss  distribution  can  be  character- 
ized by  the  most  probable  energy  loss  Q  and  the  half-width,  which  is 

P 

the  energy  width  of  the  distribution  curve  measured  at  one-half  the 
peak  intensity.   These  values  are  used  to  compare  experimental  results 
with  theoretical  predictions,  as  well  as  to  compare  with  other  work. 
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II.   THEORETICAL  CONSIDERATIONS 


The  Blunck  and  Westphal  theory  of  the  distribution  for  the  energy 
Loss  suffered  by  an  initially  monoenergetic  beam  of  electrons  in  pass- 
ing through  a  layer  of  absorbing  material  assumes  that  the  energy 
loss  Q  is  small  compared  to  the  initial  energy,  E..   Let  W(Q)dQ  be  the 
probability  of  energy  loss  between  Q  and  Q  +  dQ  and  x  equal  that 
amount  of  Q  lost  by  radiation.   Then  the  ionization  loss  equals 
(Q  -  x) .   This  leads  to  the  equation 


W(Q)dQ  =  ^  W;[(Q  -  x)Ws(x)dxdQ 


(1) 


where  W_  and  W  are  the  energy  loss  distributions  for  ionization  and 
I      s 

radiation,  respectively. 

For  the  distribution  of  energy  loss  due  to  ionization  and  excita- 
tion, W  ,  the  Landau  equation  [2],  as  modified  by  Blunck  and  Leisegang, 
is  used.   The  distribution  is  expressed  as  a  function  of  the  dimension- 
less  parameter  \: 

exp 


\-   C„Y 


W].(Q)dQ  =  $(X)d\  = 


)     n  n 
njb   +Yn 


<x  -  V 

2    2 

(bZ  +  Y„) 


dX 


(2) 


where 


^^♦-(iJ 


1.116 


(3) 


The  quantities  introduced  above  are  defined  as  follows :   c  ,  y    , 

n   n 

and  \      are  constants  used  by  Blunck  and  Leisegang  to  fit  the  Landau 
distribution  to  a  sum  of  Gaussian  functions.   Their  values  are  listed 
in  Ref.[5}  R  is  the  thickness  of  the  absorber  measured  in  cm. 
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The  quantity  a  is  a  function  of  the  atomic  number  Z,  the  atomic 
weight  A  and  the  density  p  of  the  material,  and  8  (  =  v/c)  of  the 

electrons : 

n  i  "^Zl  7  O  m0\; 

(4) 


a  = 


0.154  Z  P 


MeV 


A  8' 


cm 


The  quantity  b   is  a  correction  to  Landau  theory  given  by  Blunck 
and  Leisegang  [3~|: 

2E. 


L2   3.0 

b  =ir 


I 


m 


I  N 

m  m 


in 


I  (1  -  B") 

-  m 


(5) 


where  the  summation  is  made  over  the  m  ionization  potentials  of  the 

atomic  electrons  of  the  material,  while  N   is  the  number  of  electrons 

m 

with  ionization  potential  I  .   Blunck  and  Westphal  [ l]  and  Breuer  [6] 


used  the  approximation 


7^4/3 


b   =  2  x  10 


5  QZ 


CaR)' 


(6) 


This  approximation  is  valid  for  aluminum  and  was  used  in  making  the 
theoretical  calculations  so  as  to  permit  comparison  with  previous  ex- 
periments.  However,  the  calculations  for  copper  and  lead  were  made 
using  eq.  (5) . 

Q  is  the  average  energy  loss  due  to  ionization  (no  radiation)  for 
electrons  of  incident  energy  E.,  and  is  given  below. 

As  can  be  seen  from  eq .  (3),  A  depends  directly  on  Q.   Thus  the 
width  and  most  probable  energy  loss  corresponding  to  the  distribution 
§(\)  depend  on  Q.   Sternheimer  [7,8,9]  has  calculated  this  quantity  as 
follows  : 


Q  = 


A  t 
s 


m 


B  +  0.43  +  in   E. 

l 


:  -  c  -  afi(x1  -  log1Q£) 


MeV  (7) 
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where  t  is  the  thickness  in  gm/cm   and  the  constants  As,  B,  C,  X  ,  a  , 

and  m  are  parameters  of  the  material  and  listed  in  Ref.  9.   The  most 
s 

probable  energy  loss  due  to  ionization  and  excitation  (without  radiation) 
is  given  by : 


A  t 


A  t 
B  +  1.06  +  AriHpJ 


in 


si 


C  -  asU 


I  x   -  log   — 
\  1       10  mc/ 


MeV   (8) 


For  W  ,  the  distribution  of  energy  loss  due  to  radiation,  Blunck  and 


Westphal  give  [3,4] 


Wg(Q)dQ   =  BOftfjjp-) 


Q 


(9) 


where 


a   =   1.40  x   10 


-3   PZ' 


I  *G® + 1/» 


cm 


and  B  is  a  normalizing  factor  = 


(10) 


r(aR  +  1) 

The  distribution  of  energy  losses  according  to  Blunck  and  Westphal 
is  obtained  by  putting  Eqs .  (9)  and  (2)  into  eq .  (1)  and  performing  the 
integration.   The  results  give  an  energy  loss  probability  distribution 
for  a  single  electron  of  incident  energy  E..   For  comparison  between 
experimental  results  and  theory  this  distribution  function  has  been 
utilized,  with  corrections  to  account  for  the  finite  energy  width  of 
the  incident  electron  beam.   This  treatment  is  described  in  Section  IV. 
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III.   EXPERIMENTAL  PROCEDURE 

The  experimental  arrangement  used  was  similar  to  that  described 
by  Breuer  [6],  Miller  [121,  and  Goodwin  [131.   The  LINAC  of  the  Naval 
Postgraduate  School  was  used  to  accelerate  electrons  to  energies  from 
54  to  97  MeV.   These  electrons  were  elastically  scattered  at  90   from 
a  thin  (0.076  mm)  aluminum  scattering  foil,  thence  passing  through  the 
absorber,  and  finally  energy-analyzed  by  a  120  magnetic  spectrometer 
described  by  Kenaston,  Luke  and  Sones  [ll].   A  figure  of  this  arrange- 
ment is  given  in  Refs.  [121  and  [13], 

According  to  Miller  [12],  the  results  of  energy  loss  experiments 
with  aluminum  appeared  to  be  quite  sensitive  to  geometrical  considera- 
tions.  In  his  work  data  were  collected  with  the  absorbing  material 
positioned  approximately  22  cm  from  the  scattering  foil.   The  results 
with  this  arrangement  showed  a  general  lack  of  agreement  with  theoreti- 
cally predicted  half-widths.   Only  when  the  absorber  was  placed  approxi- 
mately 3  cm  from  the  scattering  foil  did  the  experimental  half-widths 
show  reasonable  agreement  with  theoretical  prediction  and  the  results 
of  Breuer  [6], 

Consequently,  an  experiment  was  conducted  to  determine  the  optimum 
positioning  of  the  absorber  with  respect  to  the  scattering  foil.  A 
device  was  fabricated  whereby  the  absorber  could  be  positioned  at 
various  distances  from  the  scattering  foil.   For  details  of  the  experi- 
mental arrangement  and  a  diagram  of  the  apparatus,  see  Figs.  1  and  2. 

A  portion  of  Miller's  experiment  was  repeated,  using  an  incident  energy 

2 
beam  of  53.8  MeV  and  an  aluminum  absorber  of  thickness  1.441  gm/cm  . 

The  separation  was  varied  from  3  to  20  cm  and  the  energy  distribution 
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measured.   The  results  of  these  measurements  as  well  as  the  theoretical 
prediction  are  shown  in  Fig.  3.   It  was  determined  that  the  separation 
had  no  effect  on  the  location  of  the  peak  of  the  energy  loss  curve, 
but  that  the  half-width  increased  as  the  distance  increased.   Agreement 
with  theory  as  well  as  previous  experimental  data  was  attained  for  a 
radial  separation  <  6  cm.   Thus  in  the  remainder  of  this  work  the  ab- 
sorbing layers  were  positioned  approximately  3  cm  from  the  scattering 
foil. 

Because  the  incident  electron  beam  scattered  into  the  absorbing 
material  is  not  monoenergetic ,  the  energy  distribution  of  the  electron 
beam  was  measured  both  before  and  after  passing  through  the  absorbing 
materials.   Various  thicknesses  of  each  material  were  positioned  into 
the  beam  by  a  remotely-controlled  device.   This  enabled  changing  of 
absorber  thickness  without  turning  off  the  accelerator  and  thus  poss- 
ibly altering  the  character  of  the  electron  beam.   The  data  represent 
the  number  of  electrons  detected  by  a  coincidence  counting  system  at 
the  exit  of  the  magnetic  spectrometer.   A  downstream  Secondary  Emis- 
sion Monitor  was  used  as  a  standard  for  normalization  purposes,  in 
that  each  data  point  corresponds  to  a  given  predetermined  integrated 
current.   This  corresponded  to  a  certain  number  of  electrons  passing 
through  the  absorbing  material. 
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IV.   TREATMENT  OF  DATA 

It  was  first  necessary  to  correct  all  data  taken  for  background, 
counting  loss  and  accidentals.   Background  was  measured  by  taking 
counts  with  the  scattering  foil  removed  from  the  path  of  the  incoming 
beam  and  integrating  to  the  appropriate  current.   Counting  loss 
corrections  were  applied  using  the  formula: 

n     =  n(l  +  Kn)  (11) 

corr 

where  K  is  a  characteristic  constant  of  the  counting  system  -  in  this 
case  .003  -  and  n  is  the  count  rate.   Accidental  coincidence  counts  were 
measured  by  means  of  a  63  nsec  delay  line,  and  appropriately  deducted. 
Because  of  low  background  and  the  counteracting  characteristic  of 
counting  loss  and  accidental  coincidences  these  corrections  had  neg- 
ligible effect  on  the  data. 

Since  the  theory  of  Blunck  and  Westphal  is  based  on  a  monoenergetic 
incident  beam,  the  energy  spread  of  the  incident  electrons  must  be 
taken  into  consideration.   Such  an  energy  spread  will  increase  the 
half-width,  and  because  of  the  asymmetry  of  the  energy  distribution 
will  increase  the  most  probable  energy  loss.   To  account  for  these 
effects  the  IBM  360/67  computer  of  the  Naval  Postgraduate  School  was 
used  to  unfold  the  measured  incident  energy  distribution  into  the 
theory.   The  method  is  described  in  detail  by  Miller  [121  and  Goodwin 
[l3~j.   For  a  computer  program  to  accomplish  this  unfolding,  see  Appen- 
dix C. 

In  the  comparison  of  the  data  with  theory,  the  measured  incident 
energy  distribution  was  unfolded  as  described  above,  properly  normal- 
ized to  the  experimental  data  and  plotted.   From  these  plots  the 

16 


theoretical  half-widths  and  most  probable  energy  losses  were  extract)  d. 
The  experimental  data  were  superimposed  on  the  theoretical  curves  and 
the  measurable  parameters  compared.   The  results  for  various  incident 
energies,  absorber  materials,  and  thicknesses  are  shown  in  Fig.  4 
through  23.   The  solid  lines  represent  the  theoretical  energy  distri- 
butions, while  the  crosses  indicate  experimental  data  points.   To 
avoid  clutter  on  the  figures  approximately  257o  of  the  data  points  have 
been  omitted.   A  representative  sampling  of  error  bars  is  included. 
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V.   RESULTS  AND  OBSERVATIONS 
The  theoretically  predicted  and  experimentally  measured  values  of 
half-widths  and  most  probable  energy  losses  for  the  various  incident 
energies,  absorber  materials,  and  thicknesses  are  shown  in  Tables  I, 
II,  and  III.   In  addition  to  the  experimental  data  and  the  Blunck  and 
Westphal  predictions,  the  theoretical  value  of  the  most  probable  energy 
loss  due  to  ionization,  S,  is  listed  for  comparative  purposes.   The 
data  for  an  incident  energy  of  approximately  75  MeV  for  all  three 

materials  of  various  thicknesses  are  plotted  in  Figs.  24  and  25. 

2 
For  the  thickest  lead  absorber  (t  =  2.825  gm/cm  )  at  91.37  MeV  and 

2 
the  copper  (t  =  5.726  gm/cm  )  at  53.63  MeV  a  half-width  was  not  ob- 
tained because  of  the  shape  of  the  low  energy  tail.   In  an  attempt  to 
get  a  half-width  the  energy  spectrum  of  the  electrons  was  examined 
down  to  about  5  MeV,  the  lower  limit  of  reliability  of  the  magnetic 
spectrometer.   It  was  found  that  for  thick  absorbers  secondary  effects 
(probably  due  to  shower  production)  cause  a  secondary  peak  of  low 
energy  electrons.   This  rise  began  in  the  vicinity  of  30  MeV  and  peaked 
at  about  15  MeV,  approximately  independent  of  incident  energy  and  ab- 
sorber material.   It  was  noted  that  the  amplitude  increased  with  target 
thickness,  and  in  the  case  of  the  thickest  absorbers  was  as  much  as 
ten  times  greater  than  the  primary  peak.   Measurement  of  the  presence 
of  positrons  was  made  by  reversing  the  polarity  of  the  magnetic 
spectrometer.   The  energy  spectrum  of  the  positrons  at  low  energy 
followed  that  of  the  electrons,  but  the  intensity  was  much  lower.   This 
indicated  that  the  lower  energy  tail  was  due  in  part  to  pair  production 
resulting  from  the  interaction  of  bremstrahlung  with  the  absorber 
material,  inasmuch  as  only  a  negligible  number  of  positrons  were 
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detected  when  the  absorber  was  removed.   Another  possible  contri- 
bution of  electrons  to  the  low  energy  tail  is  electron-electron 
scattering. 

From  Table  I  it  may  be  seen  that  for  lead  the  experimental  results 

show  agreement  with  theory  in  the  most  probable  energy  loss  for  all 

2 
energies  and  thicknesses  up  to  2.825  gm/cm  .   There  is  agreement  in 

2 
the  half-width  for  thicknesses  up  to  2.118  gms/cm  ,  but  at  greater 

thickness  the  experimental  half-width  is  smaller  than  that  predicted 

by  theory.   This  is  in  contrast  with  previous  works  with  materials 

of  lower  Z  where  for  thick  absorbers  the  experimental  half-widths 

are  systematically  greater  than  theory  predicts. 

2 
Table  II  indicates  that  for  a  5.574  gm/cm  aluminum  absorber  the 

most  probable  energy  loss  data  agrees  with  theory  within  the  error 

limits  of  the  experiment.   However,  at  an  incident  energy  of  both 

53.44  and  74.47  MeV  the  experimental  half-width  is  significantly 

greater  than  theory  predicts. 

For  copper,  again  agreement  between  experimental  data  and  theory 

2 
for  the  most  probable  energy  loss  was  observed  up  to  5.726  gm/cm  . 

There  appears  to  be  agreement  in  the  half-width  for  a  thickness  of 

2  ' 

4.295  gm/cm  ,  but  the  results  of  the  experiments  using  the  5.726  gm/cm'' 

absorber  are  inconclusive  owing  to  the  large  degree  of  uncertainty. 

It  would  appear  in  general  that  the  Blunck  and  Westphal  theory  for 
the  energy  loss  distribution  of  electrons  is  valid  for  predicting  the 
most  probable  energy  loss  for  aluminum,  copper,  and  lead  for  incident 
energies  in  the  range  54  to  97MeV.   The  theoretical  prediction  of  half- 
width  for  thick  absorbers  seems  best  for  copper,  but  perhaps  is  sys- 
tematically too  large  for  materials  of  higher  Z  and  too  small  for 
those  of  lower  Z. 
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In  Light  of  the  foregoing  further  work  is  justified  with  materials 
of  low  and  intermediate  atomic  number.   Furthermore,  a  closer  exami- 
nation of  the  characteristics  of  the  low  energy  tail  and  secondary 
peak  is  warranted. 


20 


APPENDIX  A 


TABLES 


TABLE  I.  Energy  Loss  Distribution  Characteristics  of  Lead 


Q  (MeV) 
P 

HW(MeV) 

E. (MeV) 

2 
t(gm/cm  ) 

B&W 

S 

Experiment 

B6W 
Folded 

Experiment 

53.85 

0.706 

0.83 

0.69 

0.81  + 

.05 

0.61 

0.60  +  .06 

1.412 

1.41 

1.43 

1.45  + 

.08 

1.21 

1.23  +  .12 

2.118 

2.65 

2.20 

2.64  + 

.10 

2.51 

2.65  +  .25 

2.825 

3.75 

2.98 

3.73  + 

.15 

7.44 

3.00  +  .30 

74.74 

0.706 

0.85 

0.69 

0.82  + 

.05 

0.83 

0.82  +  .05 

1.412 

1.80 

1.45 

1.80  + 

.10 

1.40 

1.47  +  .07 

2.118 

2.70 

2.22 

2.59  + 

.12 

2.50 

2.24  +  .19 

2.825 

3.90 

3.01 

3.80  + 

.20 

8.78 

6.03  +  1.0 

91.37 

0.706 

0.92 

0.70 

0.87  + 

.07 

0.89 

0.78  +  .08 

1.412 

1.84 

1.45 

1.77  + 

.10 

1.54 

1.52  +  .19 

2.118 

2.84 

2.23 

2.75  + 

.15 

2.85 

2.74  +  .64 

2.825 

4.04 

3.02 

3.88  + 

.20 

v- 

JL 

*  No  half-width  obtained. 
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TABLE  II.   Energy  Loss  Distribution  Characteristics  of  Aluminum 


Q  (MeV) 
P 

HW(MeV) 

E  (MeV) 

i 

t(gm/cm  ; 

B&W 

S 

Experiment 

B&W 
Folded 

Experiment 

53.44 

5.574 

9.22 

8.41 

9.24  +  .13 

3.65 

4.60  +  .30 

74.47 

5.574 

9.13 

8.43 

9.33  +  .25 

4.00 

5.7  +  1.2 

TABLE  III.   Energy  Loss  Distribution  Characteristics  of  Copper 


Q  (MeV) 
P 

HW(MeV) 

E .  (MeV) 

2 

t(gm/cm  , 

B&W 

S 

Experiment 

B&W 
Folded 

Experiment 

53.63 

4.295 
5.726 

6.28 
8.90 

5.56 
7.52 

6.23  +  .13 
8.83  +  .20 

4.30 
* 

4.25  +  1.20 
* 

76.02 

4.295 
5.726 

6.29 
9.13 

5.57 
7.54 

6.28  +  .13 
8.38  +  .25 

4.59 
9.20 

4.55  +   .50 
10.1  +  2.0 

94.07 

4.295 

6.30 

5.58 

6.40  +  .20 

4.61 

~  4 

96.66 

5.726 

9.30 

7.55 

9.2  +  .5 

15.55 

16.4  +  4.2 

*No  half-width  obtained. 


22 


APPENDIX  B 


Scattering  Foil 


FIGURE  1.   Experimental  Arrangement  for  Varying  Foil-to-Absorber 
Distance  in  Beam  Line. 
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APPENDIX  C 


Dppr-RAV    ftp    COMPUTING    ENERGY    LCSS    DISTRTBUTICN 

IMPLICIT    REAL*8     (A-H.L.Q-ZI 

OTMENSITN    E(IOO).    C7QOO),    EX(200),    CT(200) 
CHMMON/tVL/ET ,Q,QB,AR,ZtALR,BETSQ 
EXTERNAL    V* 
W«?ITF     (6,881 
pp    p^PWAT    (  *1»,11X.  'EI  STX, •BETA**2',13X, 'AP« ,11X, • ALPHA 


1    p  '  ,<=X, »OBAR«f12X, »QPI «,////) 


READ    N-NUMPER    HE    INPUT    ENERGY    BINS    C    M-NUMBER    OE    OUTPUT 
ENERGY    VALUES 

RPAO     ( 5,10)     K,M 

10  FORMAT    (2110) 
RE4"    (c,ll)    R 

11  POPMfiT    (E10.5) 

READ    7-ATO^IC    MR,    PHO-DENSTTY.    A-ATOMIC    WT ,    t    THE    PARA- 
METERS   OF    THE    MATERIAL    FROM    P*GE    514,    STFRNHEIMER    PAPER 

READ    (5,12)     7,PHO,    A,    AS,    B,    CS  ,    SMALA,    EM,    XI 

12  FORMAT     (6F10.4,/,3E10.4) 

READ  7ERO  PEAK  ENERGIES  AND  COUNTS  FROM  HISTOGRAM 

PEAO    (*.13)     <E( J) ,C7{ J) ,J=i,N) 

13  FORMAT    (F10.2,F10.0) 
Dn    ICC    1  =  1, M 

ice  ct(  n  =  c.c 

FX(1)  IS  THE  HIGHEST  OUTPUT  ENERGY  VALUE  DESIRED 

EX(  1)=SC.0C 
DLTE  IS  THE  ScPARATION  RETWEEN  ENERGY  OUTPUT  VALUES 

DLTE=.l 
COMDUTF  ALPHA*R 

ALR  =  ,CCl^*7**2*PHO*R/A*(^./3.*DLOG(18  3./7**(l./3.)  ) 
I  +  l./<3.) 

T=PHO*P 

00  2C0  J=1,N 

EI=E( J) 

COMPUTE  BETA**2 

BET SQ=1.-. 511**2 /(EI+. 51 1)**2 
COMPUTE  A*P 

AR=.15A*7*RHC*P/(A*BETS0) 
COMPUTE  P/(M*C) 

X=DSORT(EI**2-H.022*EI  )/.511 

COMPUTE    Q-BAP 

0R  =  AS*T/BETS0*(B*.43-fDL0G(EI)-BETS0-CS-SMALA*(Xl 
1    -0L0G1C( X))**EM) 

CDMPUT    Q    PROBABLE    IONI7ATION 

OPI=AS*T/BETSO*(B*1.06*-DLOG(AS*T/BETS0)-BETS0-CS 
1    -SMALA*( Xl-DLOGIO(X) )**EM) 
DO    300    1  =  1, M 
Q=EI-EX(I) 
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:X(  H-1)=FX(  i  )-rL-rr 
IFfQ.Lr.C.O)  KC    T0  30C 

UfV  ^c  THC  TNTPOMEPT  4Tf-  LIMT'C;  Oc  !NTCGR  *t  \  pm 

U=.CCCOCC1*0 
V=.2C*0 

cm.   n^r-32(0.c,i),wTFi) 

Call    nOG?2(UtVtWfF2) 

C-ML    DOr,3?(VtO,W,P3) 

p  =  F  i  +  (=2  +  F3 

CT( I )=CT( I )+F*f 7( J) 
ICC    CnNTTNjLiF 

WITr-    (f,20)    EI  »BFTSO,AR,ALRtQR»QP! 
2(\    PQPMAT    fPlS.ZfFl^^.APlS.^t/) 
2CC    rnNjTINUF 

WPT^f     (6,77) 
77    cncMU    (  »lf  ,6X,  'ENERGY1, 9Xf  •COUNTS1 1////) 

WRTtp     (^t?l)     (FX(I)fCT(I),     I=l,M) 

2i    cppM.M    (F11.2»F15.0J 

CN0 

FUNCTION    TO   CON!DUTF     the    ENERGY    LOSS    DUE    T 1    ICNIZATICN 

REAl     FUNCTION    WI  *8(FI  ,Q,GB,AR,7.  ,8ETSQ) 
IMPLICIT    R=*L*8     <A-H,L,Q-7_) 

<;Tnpn    thf    BLUNCX    f,    LFICFG6NG    CONSTANTS 

5?H*°    $(4)  /C.174,.CF8,.01°,.00  7/tG**MS  (<<  )'1.8,2.0,  3.0 
1 ,  5.  0/  v  t  A  *RDA  (4 )  /0.  0  *  3  .  0  f6  .  5  , 1 1 .  0  / 

STQRF  THP  IONIZATION  POTENTIALS  AMO  ELECTRON  LFVFLc  OF 

THF  MATERIAL 

PPAL*8    Iw(  6)  /.  0  878  9, .  01432  , .  002  <U  f .009  3  8  ,  •0000"'9, 
I     .0000  38/,NM(M  /2.  ,8.  ,18.  ,32.  ,13., A./ 

WIPN=C.C 
*R=C.C 

CALCULATE    LAVROA 

ALA M  =  ( O-OB ) /AR+OLGG ( F I  ZAP ) -1 . 1 1 6 

CALCULA^F    8**2 

DH    IOC    1*1,6 

B*IM( I)*NM(I )*CL0G(2.*FI/< I M( I ) * (1 . -BET? 0) ) ) 
100    SB  =  SB+R 

RS2=3. *$R/f 7*AP) 

nn   400   k=i,4 

PP  =  ( <  h it v-LAVftru  (K))**2) /(GAMMS (K) **?  +  R^  2) 
xw^=nFXP(-cP) 
TP( XW.IF. l.CO-60)XW=0.0 

WIN=( St K)*(GAMMS(K) /AR) /nSORT(GAMMC (  K  )  **2  +  R^  2  M*XW 
4Q0    WlON=fc!ON+WIN 

WI=WI0M 
pptupm 

END 

FUNCTION    Tn   CO^PUT^    THF    ENFRGY    LGFS    OUF    TO    Rid  AT  ION 

PFAL    FUNCTION    VF  *8  (  F I  ,  X  .  ALR  ) 

IVCLICIT    REU*8     (A-H,L,0-7) 

W*    =    AtR/DGAWMA (ALR+1. I *(X/FT I **ALF/X 

RETIJPN 

END 
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FUNCTTCN    TO   COMPUTE    TOTAL    <=NERGY    LOSS 

REAL    FUNCTION    W*8(X) 

IMPLICIT    RFAl*«     (A-H,L,Q-Z) 

CnMMON/$VL>EI,0,QB,AP,7,AL«t8ETSQ 

W=WT(F! tO-X,OB,AP,?,BETSO) *WS<EI  ,X,ALR) 

RETURN 

FNO 

SUBROUTINE    0QG32(XLfXU,FCT,Y) 

PUR°C*E 

TC   COMPUTE    INTFGRAL(FCTfX) 

USAGE 

CALL    0CG32     ( XL , XU»FCT.Y ) 

PARAMETER    FCT    REQUIRES    AN    EXTERNAL    STATEMENT 

OESCPIPTICN  OF    PARAMETERS 

XL  -  DOUBLE    PRECISION    LOWER    BOUND 

XU  -  DOUBLE    °REC!SICN    UPPER    BOUND 

PCT  -  ^UNCTICN    TO    BF    INTEGRATED 

Y  -  THE    RESULTING    DOUBLE    PRECISION   VALUE 

METHOD 

EVALUATION    IS    DONE    BY    32-PCINT    GAUSS    QUADRATURE 

DOUBLE    FPECISION    XL  ,  XU ,Y ,A ,B ,C , FCT 

A=.500*( XU+XL) 

B  =  X  U—  XI 

C=.49  863193092474078D0*B 

Y=.3  50930  500473  504  83D-2*(FCT(A-i-C)-i-FCT(A-C)) 

C=. 49  2*0 5 755772 634 17D0*B 

Y=Y+.8137197365*52  83*D-2*<FCT<A*C)*-FCTfA-C)) 

C=. 48  2381 12779375322DO*8 

Y=Y+.  1269 603265463 10300-1  *t  FCT (A+U+FCTU-C)  ) 

C=. 4674 5303 7968 86 98400 *8 

Y  =  Y+.  m 369314* 651 071 7D-1*( FCT (A+C)+FCT(A-C)  ) 
C=. 44 8160  5778 830260600*8 

Y=Y*. 21*1 794901 1113340D-1*< FCT <A*C)*FCT(A-C)  ) 
C=. 424 683806 8662 849900*8 

Y=Y*-.2  54C<3029631188088D-1*(FCT(A4-C)+FCT(A-C)  ) 
C=. 397241 897 983 °7120D0*B 

Y  =  Y+.  29 34204673 9267774D-1*(  FCT <A*CH-FCT(A-C)  ) 
C=. 366 C° 10 59370144  8400*8 

Y  =  Y+.3  29U11138  8180923D-1*(FCT<  A+C ) +FCT ( A-C  )  ) 
C=. 33 15221 334651 C76000*B 

Y=Y+. 36 172 8970544242530-1*1 FCT <  A+C  )  *-«=CT  <  A-C  >  ) 
C=. 2° 3 857 878620381 1600*B 

Y=Y*.3CC96947893  53515  3D-1*(FCT< A+C)+FCT(A-C)  ) 
C=. 25344995446611470D0*B 

Y=Y*.41655962113473378D-1*(FCT(A+C)>FCT(A-C)  ) 
C=. 210  675  6380653176700*8 

Y  =  Y+.4  3e260465  02201906O-l*(FCT(A*CH-FCT(A-C)  I 
C=. 165  934301141063  82  00*8 

Y  =  Y+. 4 55 869393*7881 942 D-l*( FCT ( A*C ) +FCT ( A-C >  ) 
C  =  .  ll<564368112606854O0*B 

Y=Y*. 4 6922 1995404022 83D-1*( FCT (A+C)>FCT(A-C)  ) 

C =.72 2 3 5980791 3 982 50- I *B 

Y=Y  +  .4  781936C03  9637430D-l*(FCT(A-i-C)*FCT  (A-C)  ) 

C=. 241538  328438691580-1*8 

Y=B*( Y+.4827C044257363900D-1*(FCT(A*C)*FCT(A-C> ) ) 

RETURN 

END 
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